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Abstract 

We present the evolution of the shell structure of nuclei in Hartree-Fock calculations using 
Skyrme's density-dependent effective nucleon-nucleon interaction. The role of the tensor part of 
the Skyrme interaction to the Hartree-Fock spin-orbit splitting in spherical spin unsaturated nuclei 
is reanalyzed. The contribution of a finite range tensor force to the spin-orbit splitting in closed 
shell nuclei is calculated. It is found that the exact matrix elements of a Gaussian and of a one-pion 
exchange tensor potential could be written as a product Skyrme's short range expression times a 
suppression factor which is almost constant for closed shell nuclei with mass number A > 48. The 
suppression factor is ~ 0.15 for the one-pion exchange potential. 
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I. INTRODUCTION 



The shell structure is a distinctive feature of nuclei and is characterized by the existence 
of magic numbers which are a consequence of the spin-orbit interaction [1, 2]. The spin-orbit 
interaction can be understood in a mean field approach which leads to a one-body potential 
containing a central part and a spin-orbit part. In spin-saturated nuclei the spin-orbit part 
stems from the spin-orbit nucleon-nucleon interaction. In spin unsaturated nuclei there are 
additional contributions coming either from the exchange part of the central two-body force 
or from the tensor force [3-5] . In view of the recent progress related to the discovery of exotic 
nuclei (neutron or proton rich) a major problem is to understand how the shell structure 
evolves from stable to unstable nuclei. Presently there is much concern about the role of 
the tensor force in the shell evolution and the structure of exotic nuclei [6-1 1] . 

In a previous work [3] we estimated the contribution of the tensor part of the Skyrme 
interaction to the Hartree-Fock spin-orbit splitting in several magic nuclei and adjusted the 
strength of the tensor force such as to obtain a good global fit. In the present paper we 
extend the previous study to exotic nuclei, most of which were unknown at that time. This 
extension sheds a new light on the previous results. 

The tensor term in the Skyrme interaction is written as a 5-function in the internucleon 
separation multiplied by momentum dependent terms (Sec. 4). The momentum dependence 
takes the finite range of the interaction into account. Contrary to the view that it plays 
a minor role because of its 5- type structure [6], this interaction has the same effect as a 
finite size interaction due to its momentum dependence. We will show that the Skyrme 
interaction provides a good mechanism for describing the evolution of the shell structure in 
exotic nuclei. 

Otsuka et al. [6] have pointed out that the nuclear tensor force has a rather long range 
and that the use of the energy density (1) may not be justified. The goal of the present 
paper is to show that expressions (1) together with (16) and (17), given in Section 2 of 
this paper, can still be used to study the contribution of finite range tensor forces, even if 
the range is that of the one-pion exchange force. Shell gaps are mainly determined by the 
spin-orbit splitting of the states with highest / in any shell and our study is restricted to 
these states. The spin-orbit splitting is less important in states with lower / because it is 
hidden by pairing effects and other forms of configuration mixing. 
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In Section 2 we derive an expression for the leading contribution of the tensor force 
to the Skyrme energy density functional using some results from the paper of Negele and 
Vautherin [12] on the density matrix expansion method. Section 3 presents some numerical 
calculations which show that the main effect of a longer range interaction is to introduce 
a suppression factor which is almost constant for all nuclei with mass number greater than 
A ~ 28. Section 4 recalls the expression of the tensor for a short range tensor interaction. 
In Sec. 5 we present results for single particle levels of Sn isotopes, N = 82 isotones and Ca 
isotopes, where the tensor force considerably improves the agreement with the experiment 
when its parameters are properly chosen. 

The conclusion is that the Skyrme energy functional with the tensor force is adequate to 
describe the evolution of shell effects. 



II. CONTRIBUTION OF A SHORT RANGE TENSOR FORCE 

The Skyrme parametrization of a short range tensor force leads to a contribution to the 
energy density 

AfT(r) = \a[J 2 n {v) + J p 2 (r)] + (3J n (r) J p (r) (1) 

where the J q (r) (q = n,p) are spin-orbit densities and a and f3 are parameters defined in 
Section 5. They represent the combined effect of the tensor plus central exchange inter- 
actions. If the radial wave functions depend only on the orbital angular momentum I and 
not on j then the spin-orbit densities are zero if both components of a spin-orbit doublet 
are filled. Then the energy density (1) brings contributions only to spin unsaturated nuclei. 
Thus AH(r) would be almost zero in 40 Ca which is a double closed shell nucleus. It would 
be large for 208 Pb which has spin unsaturated shells for both neutrons and protons. The 
energy functional (1) leads to a simple modification of the single particle spin-orbit potential 
for both protons and neutrons (see Section 5). 

The purpose of the present section is to derive the form of AiJ(r) for a short range tensor 
interaction. We focus on the contribution of the neutron-proton interaction. The starting 
point is a two-body tensor potential 

V T (r) = v T (r) f[-f 2 [ j^O^i ' ^)(°^ • r) - 6\ ■ a 2 ] (2) 
like the one arising from one-pion exchange. The effect of the isospin dependence in Eq. (2) 
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is to make (3 ~ 2a. An interaction with no isospin dependence would have (3 = 0. 

According to Negele and Vautherin [12] the expectation value of the tensor interaction 
with Hartree-Fock wave functions is 

(V T np > = - / d 3 ri d 3 r 2 Mri - r 2 ) |p n (r 1? r 2 ).&(n, r 2 ) | (3) 

The expressions for the nn and pp contributions are similar but are each multiplied by a 
factor 1/2 from the isospin dependence of (2). Negele and Vautherin give a factorization of 
the spin-density matrices for spherical nuclei in which sub-shells are either completely full 
or completely empty. It is 

p(ri, r 2 ) = i(n x r 2 )pi(ri, r 2 ) (4) 

where 

Pi(ri,r 2 ) = ±E 9 \ o RnAr^Rnj^PKcosO). (5) 

P( is the derivative of the Legendre polynomial P h 9 is the angle between the directions of 
ri and r 2 and the ± sign in Eq. (5) stands for j — I ± \-. For a short range interaction 9 ps 
in Eq. (5) and P/(cos#) 1(1 + l)/2. This is the origin of the spin-orbit splitting factor in 
Eq. (6) of Ref. [3] 

(2j + l)[j(j + l)-/(/ + l)-3/4] = ±2/(/ + l) if j = /±l/2. (6) 

If the radial wave functions R n ji(r) are the same for j — I ± 1/2 then the contribution of a 
particular /-level to pi(rx, r 2 ) vanishes if both j-components are either completely occupied 
or completely empty. 

When the interaction Vr(r) has a sufficiently short range the (V^ p ) simplifies to 

(V7) = lf dhj n {r)J p {r) j v(s)s 4 ds. (7) 

where 

J q (r) = 2rp lq (r,r) (8) 

which is the spin-orbit density in Eq. (6) of Ref. [3]. Equation (6) shows that J q (r) > 
when the lower component of a spin-orbit doublet is being filled and goes to zero when both 
components are filled. 
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III. A FINITE RANGE SUPPRESSION FACTOR 



The analysis in the present section shows that, for a tensor interaction with a range of 
the order of the one pion exchange potential and for single particle states with the largest 
I for a given A, the effect of the finite range interaction is to multiply Eq. (7) by a simple 
suppression factor which is almost the same for any nuclei with mass number greater than 
A = 28. As a consequence one should be able to parametrize the contribution of a tensor 
force to the energy density by the simple form in Eq. (1) with values of a and (5 which are 
constant for all nuclei. This means that the original short range Skyrme density dependent 
form used in [3] remains entirely valid and that finite range effects can be incorporated by 
using suitable values of a and (3 (see Section 5). 

We start with Eq. (3) and make a change of variables in the expression for (V^ p ) which 
becomes 

(Vjf p > = 8tt 2 f°° F(r)dr , (9) 
J o 

where 

F(r)=J dsrj d9 sin 3 9 v T (s)(r 1 r 2 ) 4 p ln (r 1 ,r 2 )p lp (r 1 ,r 2 ), (10) 

with r\ = r + s r /2, r 2 = r — s r /2 and 9 the angle between r x and r 2 . There are three other 
angles which have been integrated out to give a factor 8n 2 . The formula contains the factor 
\ri x r 2 | 2 = (r 2 — s 2 /4) 2 sin 2 #. The spin densities p lq are defined in Eq. (5). The squared 
distance |s| between the points ri and r 2 is 

|s| 2 = |rx - r 2 | 2 = si + 4(r 2 - S 2 /4) sin 2 t , (11) 

For states with maximum / in any shell the function F(r) has a single peak at r m near the 
maximum of the radial wave functions. The short range approximation to F(r), denoted by 
F (r), holds when the range of the interaction is much less than r rn . The important values 
of s r are much less than r m and the angle integral has contributions only from small values 
of 9. Then F(r) is replaced by 

/*oo rco 

F (r)= ds r d9 9 3 v T (s)r 8 p ln (r,r)p lp (r,r), (12) 

JO JO 

with |s| 2 = s r 2 + r 2 9 2 . The short range approximation (7) to (V^ p ) can be obtained from 
(12) by using the relation (5). 
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The numbers in Table 1 are calculated with oscillator radial wave function 



R t (r) = N t r l+1 exp 



They have a maximum at r rn = byl + 1. The radial suppression factor S(r) and the total 
suppression factor S are defined by 

Sir) S- JF{r)dr (13) 

There is a simple approximation 1(1, r /a) for the suppression factor S G (r) for a Gaussian 
interaction V(r) = Vq exp(— r 2 /a 2 ) which is given by 



1(1, r /a) - 



2/„2 \2 r 1 lV2 



1 + a 2 /6 2 



(14) 



\l(l + 1) + 4r 2 /a 2 , 

with /=(/„ + / p )/2. The form of Eq. (14) arises from the replacement of r 2 — s 2 /4 and of 
r\T2 by r 2 in Eqs. (10) and (11), an approximation which is valid when s r ~ a « r. Then 
the integral (10) separates into a product of an angle integral and an integral over s r for a 
Gaussian interaction. The first factor is independent of the radial wave functions. 

Values of 1(1, r m /a) for several closed shell nuclei are given in Table 1. The range a = 1.2 
fm is taken from [9]. There is quite a strong dependence on r for fixed / but 1(1, r /a) is almost 
constant at r = r m because 1(1 + l)a 2 /r 2 ra does not change much for all the nuclei considered. 
The table also gives values of S G (r m ) and the total suppression factor S G calculated by 
numerical integration. The approximate formula (14) for S G (r m ) is accurate to within 5%. 
The results indicate that, for a Gaussian potential and for states with maximum /, one can 
use the short range approximation (7) with a reduced interaction strength. 



A I 


r m (fm) 


I(L,r m /a) 


S G (fm) 


S G 


S Y (r m ) 




28 2 


3.08 


0.550 


0.579 


0.459 


0.197 


0.159 


48 3 


3.89 


0.515 


0.521 


0.436 


0.171 


0.147 


90 4 


4.83 


0.511 


0.507 


0.440 


0.166 


0.146 


132 4,5 


5.43 


0.520 


0.511 


0.449 


0.166 


0.147 


208 5,6 


6.34 


0.516 


0.504 


0.452 


0.164 


0.145 



Table 1: Values of the suppression factors 1(1, r m / a), S G (r m ), S G , 
S Y (r m ) and S Y for the Gaussian and one-pion exchange potentials 
for various values of A and 
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The final two columns of the table give values of S Y (r m ) and S Y for the tensor part 
of a Yukawa (one-pion exchange) potential which has a longer range and a form factor 
with radial dependence vt(x) = Vq exp(— x)(l/x + 3/x 2 + 3/x 3 ), where x = jis with ji = 
0.70 fm -1 . It has a 1/s 3 singularity at s = but this is cancelled by the sin 3 # factor in 
the integral for the matrix element. The suppression factors S Y (r m ) and S Y , calculated 
by numerical integration, are almost constant. This shows that it is reasonable to use 
the Skyrme parametrization (1) to study the contribution of the tensor force to spin-orbit 
splittings for states with maximum / even for the one-pion exchange potential. 

Our calculations show that the suppression factor S(r) is an increasing function of r 
which goes to zero as r — > and to 1 for larger r. The numbers in the table show that the 
total suppression factor S is less than the suppression factor evaluated at r rn for both the 
Gaussian and one-pion exchange potentials. 

Early calculations showed that the Yukawa one-pion exchange potential play an important 
role in describing the 208 Pb levels [13]. Note that the results given above for the Gaussian 
potential shape are important in view of the fact that such interactions are used in shell 
model calculations (cf. Otsuka et al. [9]). The main difference between the values for 
S (r m ) and S Y \r m ) in Table 1 is due to the range of the Gaussian interaction. The values 
become very similar for all the nuclei in the table if the range of the Gaussian interaction is 
increased from 1.2 fm to 2.1 fm. 

IV. THE TENSOR PART OF THE SKYRME INTERACTION 

The parameters of the Skyrme interaction were originally determined in Hartree-Fock 
calculations to reproduce the total binding energies and charge radii of closed-shell nuclei 
[4]. Further extensive calculations were made later [5]. Several improved parameter sets 
were found. They differ mainly through the single particle spectra. In the present paper 
as in our previous work, we shall use the parameter set SIII which gives good overall single 
particle spectra. In Ref. [3] a tensor force was added and a range of its strength was found 
such as to maintain a good quality of the single particle spectra of 48 Ca, 56 Ni, 90 Zr and 208 Pb. 

As in Ref. [3], in the configuration space the tensor interaction has the following form 

V T = ^T{[{6\ ■ k'){a 2 ■ 9) - ^'Vi ' fiMfi - f 2 ) 
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1 

+S(f[ - r^)[((fi • k)(a 2 ■ k) - -k (&[ ■ <T 2 )]8(f[ - f[)} 

+U{(a 1 ■ k>)5(f[ - r^a, ■ k) - ^(a, ■ a 2 )[k> ■ 5(n - r 2 )k]}. (15) 

The parameters T and U measure the strength of the tensor force in even and odd states of 
relative motion. 



V. RESULTS 

The analysis in Sections 2 and 3 show that the simple form (1) is a good approximation 
to the contribution of the tensor forces to the energy density. Values of a and (3 can be 
taken to be constant for states with maximum / in nuclei with A > 48 even for forces with 
a range of the one pion exchange potential. 

Both the central exchange and the tensor interactions give contributions to the spin- 
orbit single particle potential to be added to the spin-orbit interaction. The additional 
contribution are [3] 

AW n = (aJ n + pj p )e-s (16) 

AW p = (aJ p + (3J n )£-s (17) 

with a = ax + ol c and (3 = 0t + Pc- F° r the Skyrme SIII interaction the parameters of the 
central exchange part are [5] 

a c = hh - t 2 ) = 61.25 MeV fm 5 , C = , (18) 
8 

where t\ and t 2 are two of the Skyrme interaction parameters. In terms of the tensor 
parameters T and U one has 

a T = ^U, fr = l-(T + U). (19) 

Equations (16) and (17) imply that the mechanism invoked by Otsuka et al. is intrinsic 
to the Skyrme energy density formalism. These equations show that the filling of proton 
(neutron) levels influences the spin-orbit splitting of neutron (proton) levels whenever f3 ^ 0. 
The normal spin-orbit single particle potential is 

V so = W,\ d f + ^)i-s with ^<0. (20) 
r ar ar ar 



When (3 is positive the neutron (proton) spin-orbit splitting is reduced as protons (neutrons) 
fill a j — I + 1/2 level because J p(n ) > 0. This effect is clearly seen in Fig. 4 of Otsuka et al. 
[7]- 

In Ref. [3] we searched for sets of parameters a and (3 which simultaneously fit absolute 
values of single particle levels in the closed shell nuclei 48 Ca, 56 Ni, 48 Zr and 208 Pb. We 
found that the common optimal values were located in a right angled triangle with sides 
a = - 80 MeV fm 5 , (3 = 80 MeV fm 5 and hypotenuse a + (3 = 0. Here we relax these 
constraints and try to analyze single particle energies some nuclei far from the stability line. 
The experimental data did not exist in 1977 when we discussed the global fit for closed shell 
nuclei [3]. Our present choice of parameters is guided by the recent results of Ref. [10] on the 
Z = 50 isotopes and N = 82 isotones which were analyzed in a HF + BCS approach based 
on the Skyrme interaction SLy5 [14] with refitted values of T and U plus a pairing force. To 
see whether or not one can obtain the correct trend in the evolution of single particle levels 
we look at energy differences between them. These differences can give a clear indication of 
the formation of closed shells from the size of the gaps. Absolute values of single particle 
energies depend not only on the tensor, but also on other parts of the Skyrme interaction. 
Here we are not concerned with making the best fits to absolute energies. 

In the present paper we still use the SIII version of the Skyrme interaction [5] for 
comparison with the previous work. We maintain the conditions a < and (3 > but take 
values outside the triangle found before which are not inconsistent with the previous findings 
[3]. We show that values a T = - 180 MeV fm 5 and (3 T = 120 MeV fm 5 , or equivalently a = 
- 118.75 MeV fm 5 and (3 = 120 MeV fm 5 , give a reasonably good fit to Z = 50 isotopes and 
N = 82 isotones. These values are similar to the ones fitted by Brown et al. [8] in a recent 
paper. For a more general orientation we also discuss the role of this parametrization on Ca 
isotopes. 

We conclude this section with some remarks on 208 Pb and 90 Zr. The proton hn/ 2 and 
neutron i 13 / 2 in 208 Pb are filled and J p and J n are both positive with comparable magnitudes. 
Because a xs —(3 we have AW n xs AW P w and the tensor forces hardly change the spin- 
orbit splitting. The situation is different for 90 Zr. There J p = and the effect of the tensor 
forces is to increase the spin-orbit splitting for neutrons and reduce it for protons. The shell 
gaps for protons and neutrons are both increased significantly and the stability of the double 
closed shell at 90 Zr is enhanced. 
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60 65 70 75 80 85 

N 

FIG. 1: The proton single particle energy difference between l/in/2 and I57/2 in Sn isotopes (Z = 
50, N = 64-82) calculated without and with tensor force a = - 118.75 MeV fm 5 , f3 = 120 MeV fm 5 . 
Data points are from Ref. [15]. Solid dots give information from transfer reactions. Open circles 
are obtained from methods less sensitive to the single particle nature. The parentheses indicate 
less certain or indirect assignments. 



A. Sn isotopes 

Fig. 1 shows the HF results for the proton single particle energy difference between 
1^11/2 an d 1^7/2 m Sn isotopes (Z = 50, N = 64-82) with and without tensor force. One 
can see that the effect of the tensor force is indeed important. The experimental pattern 
is satisfactorily reproduced with this simple approach. In the more sophisticated HF+BCS 
calculations of Ref. [10] the theoretical results beyond 126 Sn are better. However in that 
region the experimental situation is less certain because the corresponding values have been 
assigned using methods which are less sensitive to the single particle nature of the levels [15]. 
For the double magic nucleus 132 Sn the effect of the tensor force and of the central exchange 
part cancels out because J p ~ J n . For isotopes with Z = 56-62 the comparison with the 
experiment is not possible because the l/in/2 level becomes unbound in these calculations. 



10 





i i 

N = 82 isotones 


i 




i 




1 


3 


no tensor - ~ - 





- - ~ 








2 














1 


(°) 















• 


• 


• 








1 


with tensor 
i i 


i 




1 




1 




50 55 


60 




65 




70 



z 

FIG. 2: The neutron single particle energy difference between H13/2 and I/19/2 in N = 82 isotones 
calculated with and without tensor force. Data points are from Ref. [15]. Solid dots give informa- 
tion from transfer reactions. Open circles are obtained from methods less sensitive to the single 
particle nature. The parentheses indicate less certain or indirect assignments. 

B. N = 82 isotones 

In Fig. 2 we present neutron single particle energy differences between H13/2 and I/19/2 
in N = 82 isotones calculated with and without tensor force and compare them with data 
from Ref. [15]. Again the role of the tensor force is considerable, bringing down the energy 
difference e(lzi3/2) — e(lhg/2) close to the best known experimental values. For Z < 50 and 
for Z > 70 the I/19/2 level becomes unbound both with and without tensor force. 

C. Ca isotopes 

As the Skyrme interaction SIII was fitted to closed shell nuclei it has the peculiarity 
that the central exchange interaction produces some undesirable effects in the middle of a 
shell. The predicted single particle levels have the wrong order when compared with the 
experimental levels and wrong levels are occupied [5] . This happens in the absence of the 
tensor interaction, but when a tensor interaction with adequate parameters is added the 
problem is solved. In particular the parameters a = - 118.75 MeV fm 5 and (3 = 120 MeV 
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FIG. 3: The proton single particle energies in Ca isotopes relative to \d2>/2 level, calculated with 
tensor force {nlj + T): a = -118.75 MeV fm 5 , j3 = 120 MeV fm 5 and without tensor force (nlj): a 
= 61.25 MeV fm 5 , (5 = 0. Data points are from Ref. [16]: solid dots for l<i5/2 and stars for 2sl/2. 

fm 5 remove this anomaly in 50 Ca. The reason is a considerable increase of the spin-orbit in 
the 1/ shell which shifts the I/7/2 above the 2p levels. However the anomaly persists for a 
= and P = 80 MeV fm 5 located at the edge of the above mentioned triangle. 

In addition, from Fig. 3 one can see that the effect of the tensor interaction is important 
and improves the spin-orbit splitting in the Id shell. In the 2si/2 shell the trend is correct 
but the theoretical results are above the experimental points, with or without tensor. The 
pattern is quite similar to that obtained in Ref. [6] in a shell model approach. 

VI. CONCLUSIONS 

The short range approximation for the contribution of a tensor force to the spin-orbit 
splitting in nuclei was studied in sections 2 and 3 for both a Gaussian and a Yukawa (one- 
pion exchange) interactions for states with maximum / in any shell for several nuclei with 
mass number between A = 28 and A = 208. It was shown that the exact matrix elements of 
the one-pion exchange tensor potential could be expressed as a product of the short range 
expression (7) and a suppression factor S Y ~ 0.147 which is almost constant for nuclei with 
mass number A > 48. It is only slightly larger, i. e. S Y ~ 0.16 for nuclei near 28 Si. Thus 
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the short range formulae (1), (16) and (17) with constant a and (3 should give qualitatively- 
good results for a Yukawa one-pion exchange potential. 

We have made a new fit to the parameters a and (5 in the parametrization (16) and 
(17) of the tensor contribution to the spin-orbit coupling using data on Z = 82 isotopes 
and N = 82 isotones. The tensor force makes a dramatic difference to the single particle 
energy difference between the hn/ 2 and g7/ 2 single particle levels. A similar situation holds 
for the energy difference between the ii 3 / 2 and h 9 / 2 single particle levels in N = 82 isotones. 
In both cases the calculation with the addition of the tensor force give a good description 
of the experimental data. The case with Ca isotopes is similar to 90 Zr. The tensor force 
reduces the spin-orbit splitting for p rotons and increases it for neutrons. This brings the 
order of single particle into a better agreement with experiment. 

The mechanism observed by Otsuka and collaborators [6, 7, 9] that the filling of neutron 
levels influences the proton spin-orbit splitting and vice-versa is intrinsic to the Skyrme 
energy density approach and is very simple in that theory. In addition, with the Skyrme 
density formalism, one can easily study the combined contribution of the central exchange 
and tensor NN interactions to the spin-orbit potential. 
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